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Abstract 

We consider double-inflationary models with two noninteracting scalar fields, 
a light scalar field (pi with potential ^mf(j)f and a heavy scalar field cj)h with 
potential y^(f>^ with n = 2,4. CDM with the initial spectrum of adiabatic 
perturbations produced in these models is compared with observations. These 
models contain two more free parameters than the standard CDM model with 
an initial scale-invariant spectrum. We normalize our spectra to COBE DMR 
and compare the predictions with observations on the biasing factor, large- 
scale peculiar velocities, quasar and galaxy formation and the Stromlo-APM 
counts- in-cells analysis. The model with n = 4 is excluded by the data while 
for the n = 2 model, taking cosmic variance into account, a small window of 
parameters compatible with observations is found. 

PACS Numbers: 
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INTRODUCTION 



Inflationary models |I| can solve some of the outstanding problems in cosmology. Also, 
as emphasized already some time ago p[, for a given model it is possible to calculate the 
various spectra of fluctuations produced during the inflationary phase and to put the pre- 
dictions to test using observational data. In the simplest inflationary models, the Fourier 
components of the gravitational potential are Gaussian stochastic quantities with an ap- 
proximately flat (Harrison-Zel'dovich) r.m.s. spectrum For more complicated spectra, 
calculation of the spectrum can be implemented either analytically or numerically, (for the 
latter possibility, see for ex. 0), and the spectrum can be obtained with very high accuracy. 
The (approximately) flat (Harrison-Zel'dovich) spectrum, together with the assumptions of 
the standard CDM (Cold Dark Matter) model was put to test using A^-body simulations 
and a very good agreement was found with the observed galaxy-galaxy correlation function 
on scales (0.5 - 10)/i-iMpc when h ^ 0.5 JT], see also ^ (h = Hq/IOO km/s/Mpc). Obser- 
vations seem to imply more power on scales greater than 20/i~^Mpc, strong evidence coming 
from the APM galaxy survey and observed peculiar velocities and from large-angle ^ 
fluctuations [ITOHT^. There are several ways to solve this problem, either abandoning the 



assumptions of CDM, one may consider a mixture of hot and cold dark matter |jT3| or a 
nonvanishing cosmological constant, or trying CDM with an initial perturbation spectrum 
having more power on large scales. Following the second possibility, an attractive solution is 
to consider so called "tilted" spectra, with spectral index n < 1. Such a spectrum will occur 



in power-law and in extended inflation. Recent studies however |14,15| have shown that no 



value of n seems able to reconcile the CDM model with all the observations (in accordance 



with the earlier "pre-COBE" discussion in [|T^,0), although the value n ~ 0.7 — 0.8 comes 
closest to it. 

So, if we want to reconcile the CDM model with observations without introducing neu- 
trinos with a restmass of a few eV, we have to consider models belonging to the next level 
of complexity, i.e. having at least one more additional parameter characterizing the initial 
perturbation spectrum. Inflationary models with flat spectra need one free parameter to 
specify the amplitude of the spectrum of density perturbations, for example the coupling 
constant of the inflaton in the simplest versions of inflation, while those generating a tilted 
spectrum need one more free parameter, for example the Brans-Dicke parameter u in ex- 
tended inflation models. We will consider here double inflationary models |ll7|-p5[]. In the 



speciflc models considered here, the inflationary stage is driven by scalar flelds without mu- 
tual interaction potential. As a general rule, such models produce a spectrum of adiabatic 



perturbations having more power on scales larger than some characteristic scale p2[. Three 



free parameters are now needed: one for the height and the form of the "step" , one for its 
location and flnally one for the overall normalization. We study here two double inflationary 
models consisting of two noninteracting scalar flelds, a light scalar fleld (pi with potential 
^mfcpj and a heavy scalar field (ph with potential ^0^. The spectra obtained numerically are 
then compared with observational data. Recently, an analogous study was performed for 
another double-inflationary model |^ and we will present the results in a way which makes 
the comparison easier. In section I we give some basic results for the models considered here 
and specify their free parameters. In section II, we give the various observational tests to 
which our models are submitted for different values of their parameter. In section III flnally, 



2 



we give a brief discussion of our results. 



I. THE MODEL AND ITS FLUCTUATION SPECTRUM 

Let us start with a short description of the homogeneous background. We consider the 
following Lagrangian density describing matter and gravity 

L = + \Mt - + \{Mr - 0?) (1) 

where = 0, .., 3, c = ^ = 1 and the Landau-Lifshitz sign conventions are used. When 
n = 2, A = m|, whereas for n = 4, A is a dimensionless coupling constant. The space-time 
metric has the form 

ds^ = dt^ - a\t)6ijdx'dx^ , i,j = 1, 2, 3. (2) 

Spatial curvature may always be neglected because it becomes vanishingly small during the 
first period of inflation driven by the heavy scalar field. The homogeneous background is 
treated classically, it is determined by the scale factor a{t) and the two scalar fields 4'h,4'i- 
Their equations of motion are given by 

d = aH, = ^i'Pl + 0? + 2-rH + 

3 n 

4>h + 3H<Ph + X<Ph~' = 0, 4>i + ^H^i + m^<Pi = 0, 

where a dot denotes a derivative with respect to t. We have the useful equation 

H = -AnGi<Pl + <pf) (3) 

which shows that H always decreases with time in this model. 

Let us consider now the inhomogeneous perturbations. We consider a perturbed FRW 
background whose metric, in the longitudinal gauge, reads 

ds'^ = (1 + 2<^)dt^ - a2(t)(l - 2m)5ijdx'dx^ (4) 

(in Bardeen's notations [^, $ = ^a-,^ = ~^h)- We get from the perturbed Einstein 
equations (exp(ikr) spatial dependence is assumed and the Fourier transform convention is 

$ = ^, 

54>i + + + m])5cPi = 40;$ - 2m20;$ (5) 

We see that, contrary to the case when only one scalar field is involved, when we have 
more than one scalar field, the dynamics of the perturbed system cannot be described by 
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just one equation for the master quantity $ or else for the gauge-invariant quantity 
C, = 5(j) + in terms of which the action for the fluctuations can be written [29 . 

Analogously to one can now compute the spectrum of growing adiabatic perturbations. 
These perturbations arise from the vacuum fluctuations of the scalar fields iph and 0;. The 
fluctuations are Gaussian and the power spectrum of the gravitational potential, 

defined through = $^(A;)5(k — k') characterizes them completely. For scales crossing 

the Hubble radius when both scalar fields are in the slow rolling regime, the spectrum of 
growing adiabatic perturbations, when those scales are outside the Hubble radius during the 

2 

matter-dominated stage (assuming a{t) oc ts at the present time), is given by 



ki^(k) = l^T,GH\H^f + 02 (6) 




4 2 

5 V n 



= ( )4^1n4-I k<^kf (7) 

on k 

where the r.h.s. has to be taken att = tk, the time at which a perturbation with wavenumber 
k comes outside the Hubble radius during inflation, k = a{tk)H{tk). The wavenumber kf 
corresponds to the characteristic scale appearing in the spectrum and it is close to the scale 
crossing the Hubble radius near the end of the first inflation. Hence, there exists a very 
broad interval of scales for which the dynamics of inflation at the time of the first Hubble 
radius crossing is determined by the field (ph while the main contribution to 0(/c) is made 
by 0/. We see from (|^ that the upper part of the spectrum, corresponding to small k's or 
very large scales, is not flat but has a logarithmic dependence cx In* This gives, as we 
will see, a crucial difference between the n = 2 and n = 4 spectra. The quantity s{t) is the 
number of e-folds from time t till the end of the second inflation and it is given by 

s ^ 4vrG(^ + f ) (8) 
n 2 

In order to have a sufficiently long second inflationary stage that will put the characteristic 
length scale of the spectrum on a scale in agreement with observations, we have that 0; ~ 
3Mp near the end of the first inflation so s{tf ) ~ 60. A very little change in this initial 
value will be enough to shift the spectrum in fc-space while leaving the form of the spectrum 

practically unaltered. We define the parameter p = ^'^^^ i gives the order of magnitude 
of the ratio of the energy of the heavy scalar field to the energy of the light scalar field near 
the end of the first inflation; it specifies the form of the spectrum, namely its "step" with 
more power on large scales, and the width of the transition region in fc-space. We will now 
estimate A^, the height of the "step" in the spectrum between a scale which is on the upper 
plateau (but of course still inside the cosmological horizon) and a scale at the beginning 
of the lower plateau. For scales crossing the Hubble radius at the beginning of the second 
inflationary stage (neglecting possible small oscillations), we have the standard result of 
inflation driven by one scalar field in slow-rolling regime 



24:TiGmf 



k^^{k) = ^ ^so (9) 
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-VQTi^G^micPl (10) 





where Sq = 27iG(j)l. We therefore get for 



Afe ^ ^ 7 (11) 



mi 

H k 



yX4^^J^ (12) 



For n = 2, this is 



mi n ^ y/so 



Afc ~ 0.13j9ln5 -1 k<t:kf (13) 
~ O.33p0;,(4), (14) 



whereas n = 4 yields 



k 

Afc ~ 0.073J9 In A; < A;/ (15) 

rh 

~O.24p0V(tfc). (16) 

where (l)h = 4r- We will investigate models with 3 < n < 16 for n = A and 6 < p < 28 for n = 
2. How we choose the scales of our spectra is very important when we compare them with the 
observations. For this purpose, we need a precise definition and we adopt here the convention 



of |2^, and define it with the help of fcb, the scale where the extrapolated upper part 
intersects the lower plateau. One shows numerically that the evolution of the background 
also during the transition between the two main inflationary stages is inflationary, in the 



sense that a > 0, for p < 25 when n = 2 |Q and for p < 50 when n = 4. Also, the initial 



perturbation spectrum will have no oscillations for p < 15 when n = 2 and also for all p's 
considered here when n = 4. A last comment concerns the power spectrum P{k) defined 
through (SkSl,) = P{k)6{k — k'). Linear perturbations grow at different rates depending on 
the relation between their wavelengths, the Jeans length and the Hubble radius and this is 
specified by the transfer function T{k) [more accurately, one should write T{k,to)] 

P{kM = l^^\k)T\k). (17) 

where by definition, T{k 0) = 1. T{k) is computed numerically making assumptions 
about the matter content of the universe, and depends on parameters like VLq and h. We 
use here the transfer function for the standard CDM model given by 

T(g) = Mi_L|:M[i + 3.89g + (le.lg)^ + (5.46^)=^ + (6.72g)^]-3 (18) 

where q = Q^f^2^.jp^-i , = I, h = 0.5. We will assume tacitly in all the formulas that 
the density parameter Qq = 1, in accordance with inflation. The power spectra P{k) for 
different values of the parameters are displayed for the n = 2, resp. n = A model in Fig.l, 
resp. Fig. 2. 
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II. CONFRONTATION WITH OBSERVATIONS 



A. Normalization of the spectrum to COBE DMR 

Let us start with the normahzation of the spectrum. Since the discovery of large angular 
scale fluctuations in the Cosmic Microwave Background Radiation by COBE DMR |]Tn|, one 
can normalize the spectrum of fluctuations using these COBE DMR observations. These 
observations probe the spectrum of fluctuations on very large scales, from several hundreds 
Mpc up to the cosmological horizon and allow for a normalization based on first principles. 
More precisely, one has 

2/ -I- 1 

4(10°) = Y. -^{\aim?) exp[-/(/ + l)el] ^ ^ exp[-/(/ + l)^^], (19) 

l>2 ^'^ l>2 

where = -4:259 fwhm = 4.25° is the Gaussian angle corresponding to the antenna beam 
and additional smearing of the raw data. It is to be noted that due to the exponential factor, 
only the multipoles up to / < 20 contribute significantly to the observed anisotropy. The 



r.m.s. coefficients ^{\aim\'^), with aim defined by 

AT 



O-lmYlm (20) 



Lm 



are actually independent of m. Based on the latest results of COBE DMR |11,12|, we take 
a^{10^) = (1.25 ± 0.2) X 10~^, with error bars at the la-level. This allows us to normalize 
the fluctuation spectrum and constitutes a great step forward. For the CMBR anisotropy, 
on large angular scales the dominant effect is the Sachs- Wolfe effect and, for adiabatic 
perturbations, we have the fundamental relation 

iM') = ^ / dkk~'Pik)jfikrrec), (21) 
Ztc Jq 

where ji is a spherical Bessel function and Tree is the comoving distance between us and the 
surface of recombination, we have in very good approximation r^ec = "J^- the models 
we are considering here, normalization of the power spectrum P{k), which is obtained for 
given p and location of the "step" , through eq. (^I]) will fix the value of the remaining free 
parameter A or equivalently rrih. The parameter rrih will have the order of magnitude given 
thereafter when we vary p and ki, earlier [^. As a result, the energy density at the beginning 
of the second inflation is also fixed: it is ~ ^ x 10~^^Mp for n = 2. For the parameter A, we 

get the following result, when we vary the parameters: -\/A ~ 2 x 10^^ and the corresponding 
energy density is then ~ ^ x 10~^^Mp for n = 4. Another remark concerns the contribution 
to the CMBR anisotropy on large scales (/ < 40) which comes from gravitational waves 
(tensor metric perturbations). One can show that this contribution is equal for both models 

and rather small, namely ^ {\aim\'^)tot = y^(l + |)(|a;m|2)^p 1.05^(|a;mp)AP where the 
subscript AP refers to that part of the fluctuations due to adiabatic perturbations, i.e., the 
quantity calculated in (^) . Although this effect is not large, it is important here for the n = 2 
model, which for some choices of the parameters will turn out to be in marginal agreement 
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with observations. Also, one has to take into account "cosmic variance" which is connected 
to the fact that we try to estimate r.m.s. values of physical quantities in the universe from 
a limited sample. Therefore, independently on how good the COBE measurements may be, 
the quantities Y^-m ^[^^"'[2^ obey a x^-distribution for 2Z + 1 d.o.f. 

In connection with recent interest in relations between ^ = l,"'"" li*^^ (for 1 < / < 40) 

and the slopes of power spectra of adiabatic perturbations, = 1 + , and of 

gravitational waves, tit = ^^^^^^5^^^; it is interesting to note that the relation ~ — 1 
(proposed, e.g., in P^J53| ) is valid in our models for k kf (and not too small) though 



it is not valid for k > kf and also not valid for single chaotic inflation. On the other 
hand, the relation ^ ^ 6.2nT which is valid for single slow-rolling inflation (see e.g. ||3^ for 



chaotic inflation) is strongly violated if khor kf (this issue will be addressed in a separate 
publication 



B. Large-scale peculiar velocities 

A lot of information can be gained from the observation of peculiar velocities, velocities in 
addition to the Hubble flow. As all matter contributes gravitationally, the peculiar velocities 
sample all the mass and not just the galaxies. Hence, knowing the peculiar velocity field 
would give us an information on the primordial spectrum of the same interest as the CMBR 
anisotropy. It should however be stressed that large-scale peculiar velocities have rather 
large uncertainties. In the linear regime, gravitational instability produces a velocity field 
that is irrotational at sufficiently late times, the velocity field v then derives from a velocity 
potential \l/ with v = — V\E'. Measurements of redshifts z and of galaxy distances r (actually 
of Hor) provide the radial component of the peculiar velocity field:fr- = cz — Hqt, for a galaxy 
at small redshift. For a potential fiow, the radial velocity field when integrated along radial 
paths gives the velocity potential out of which the other velocity components can be derived. 
The difficulty with this method is to construct a smooth radial velocity field and to eliminate 
the statistical uncertainties in Vr- This is done with interpolation and smoothing of the raw 
data and we finally have the following equation 

(^')^ = ^2 / dkP{k)W^H{kr) eM-kRs) (22) 

ZTT^ JO 

where WTH{kR) stands for the Fourier transform, up to a constant, of the "Top Hat" window 
function 

W{kR) = (sin kR - kR cos kR). (23) 

[kRy 

We will compare our predictions with the Potent data. In these observations the raw data 
are smoothed with a Gaussian smoothing radius Rs = 12/i~^Mpc while spheres of radii 
R = 40/;.~^Mpc and R = QOh~^Mpc were considered. The data are ^40 = 405 ±60 km/s and 
^60 = 340 ± 50 km/s ||36|, error bars at the la level. Analogously to what was found in [p6| , 
this test is crucial for our models too. 

Velocities are generally too low and all the models are excluded if one doesn't invoke 
cosmic variance for the measured peculiar velocities. For the n = 2 models, the best velocities 
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are systematically obtained for ~ 6/i~^Mpc. This is also the case for n = 4, except when 
p < 6, but the velocities do not grow significantly for smaller kb- Also, it turns out that the 
n = 4 model can be excluded: for p < 5 one gets the best velocities, but other tests exclude 
these models while for p > 5, the velocities become too small. 



C. The biasing factor b 

Before the COBE DMR observations of the CMBR anisotropy, one way to normalize 
the fluctuation spectrum was through the quantity cxg = ((^)^)_R=8/i-iA/pc- This quantity 
measures the variance of the total mass fluctuations in a sphere of radius R = 8h~^Mpc. The 
reason for considering spheres of radii R = 8h~^Mpc is that for bright galaxies cxg is equal to 
P^. However, one doesn't expect the total matter to be as clustered as bright galaxies 



one 



and one tries therefore the simplest assumption, namely that there is a scale independent 
bias, given by the biasing factor b 



a 



R,9 



b'cr'n = b'^ (24) 



where the subscript g refers to galaxies while ^ is the two-point correlation function. ct|. can 
be computed from the power spectrum P{k) 

1 

^1 = ^2 1 dkk'W^H{kR)P{k). (25) 

Early numerical simulations of CDM models were able to reproduce the correct two- 
point correlation function ^ which is certainly a great success for CDM. It required however, 
if one imposes Qq = 1, that h is as low as /i ~ 0.25 . However for h ~ 0.5 one gets the right 
amplitude for ^ assuming that 6^2. Another reason for considering a bias comes from the 
observed r.m.s. peculiar velocities between galaxy pairs. If f2o = 1 and there is no bias then 
numerical simulations give velocities ~ lOOOkm/s, as expected also on theoretical grounds. 



much larger than the observed ones ~ 300±50km/s (see however for a possible velocity 
bias). Here also, an f2o = 1 model will be compatible with the observed peculiar velocity 
fields if we introduce a biasing parameter b ~ 2.5. Finally, b = 2.0 ± 0.2 is needed in order 
to get the correct amount of clusters of galaxies p9| . 

For n = 4, we get an unacceptably low b {b < 1.5) for p < 5, for p = 5, 6 ~ 1.5. For 
n = 2, we get too low 6 for p < 8 which are just the values that produce the best velocities. 



D. Formation of galaxies and quasars 

The very existence of compact objects observed at high redshifts constrains any proposed 
model for the formation of galaxies as these objects must already have formed at these high 
z. Quasars have now been observed up to redshifts near z = 5. They are believed to be 
powered by massive black holes located at the center of galaxies. From luminosity bounds 
one can estimate the mass of the black holes (M ^ lO^M©), while for the host galaxy 
the estimates give M ^ 10^^ — IO^^Mq. Although the formation of gravitationally bound 
objects is a complicated non-linear process one can, using rather simple assumptions, make 



a connection with the linear theory ||8| j40|j4l[| . The fraction -F(> M) of bound objects with 
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mass greater than some given mass M can be expressed as a function of an |^^, where R 
is the radius of the sphere containing an amount of mass M today, and is calculated 
assuming a linear evolution (in this subsection we will adopt the notation a{M): 

F{> M) = erfc(-^) (26) 
v2 

where erfc is the complementary error function, u = 5c{1 + z)(t~^{M). The spherical collapse 
model gives for the collapse threshold value 6c = 1.686, the value used here, though estimates 
from numerical simulations suggest other acceptable values (1.33 < 5c < 2). The most recent 
estimates of the mass fraction F{10^^Mq) in host galaxies of quasars at z = 4 yields 



the following lower bound |2^, using ( p6D : 

(t(IO^^Mo) ^ 2.2 ± 0.5 (27) 



When comparing our predictions with the data, we should keep in mind that equation 
is only a lower bound as quasars do not necessarily form in all potential host galaxies and 
the real number of quasars at z = 4 might be larger than the observed one. Also important 
is the fact that many large galaxies seem to have formed already at z = 1. From this 
observation we get the lower bound 

ailO^^M^) ^ 2.0 ± 0.4 (28) 



Approximately the same estimate follows from Ref. |43 



For n = 2, p> 25, considering the best velocities obtained for |^ ^ 6h~^Mpc, the values 
obtained for this test are too low. For these models, it is interesting to point out that if the 
overall normalization goes up, hence improving these numbers, the biasing factor of these 
models will become too low so that these models must be rejected. 



E. Counts-in-cells analysis 

We finally compare our models with the counts-in-cells analysis of large- scale clustering 
of the Stromlo-APM redshift survey [44]. Values for the counts-in-cells variance af, where 



/ is the cell size expressed in h~^Mpc, obtained with our spectra normalized according to 
(Tg = 1, corresponding to optical galaxies in redshift space, are compared with the Stromlo- 
APM data. In order to decide whether our model fits the data well, we apply a analysis. 
Considering the 9 data points (for 9 different cell sizes) as independent and the error bars 



quoted in |Q as 2a ones, while we test here a theory with 2 parameters p and kb (we still 
have the possibility to change the normalization, this is just changing 6 ), we have a 
distribution with 7 d.o.f. The variance af can be written as 

af = — / dkeW^{kl)P{k). (29) 

where is the analogous of W^^^, but now for a cell of size / 

W^{kl) = 8 C dx C dy dz{l - x){l - y){l - z f^^^^'''' r = \x\. (30) 

JO Jo Jo rCLV 
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A < 7 will be considered good while > 18 will be considered bad. 

The n = 4 model gives very bad numbers, > 30 for p < 5, > 20 for p < 8, for 
p = 8 the test is still not too good, > H- All the models with n = 2, p > 10 will yield 
very good results, 2 < < 3 for |^ ^ (6 — 10)/i"^Mpc. 

III. DISCUSSION AND CONCLUSION 

The obtained large scale peculiar velocities are low for all values of p and fcf,. A possible 
solution to this problem is to assume that, due to "cosmic variance", the velocities observed 
around us are higher than the real average. If we assume that the peculiar velocities are 
measured from just one independent volume, v"^ (actually ffg or vl^) itself obeys a x^- 
distribution (with one d.o.f. and variance (T^(f ^) = 2), for example there is a 20% probability 
to have f ^ > 1.6(f^) and a 30% probability to have v"^ > l.l{v'^). But one has to be cautious, 
when invoking cosmic variance since a x^-distibution with one d.o.f. for v"^ might be too 
crude. 

When n = 4, we get the highest, though still rather low, velocities for the lowest values of 
p. But for p < 5, we get unacceptably low b and unacceptably high values for ^(IO^^Mq) and 
a{10^'^ Mq). Comparison with the counts-in-cells analysis gives very bad results for p < 8. 
For p = 6, 8 we get acceptable ^(IO^^Mq) and ^(lO^^M©), however the velocities become 
unacceptably low and this situation becomes only worse with growing p. For example p = lA 
gives unacceptably high h and unacceptably low ^(lO^^M©), ailQ^'^Mo) and v. Hence, the 
n = 4 model is very unlikely and can be confidently excluded. 

An improved situation is obtained when n = 2. We obtain a window of allowed param- 
eters for 1^ ~ (6 - 10)/i-iMpc and 10 < j9 < 15 (see fig.3). However, also for this window, 
the velocities are too low compared to the POTENT data and one has to invoke "cosmic 
variance" and higher observed velocities around us than the real average one. In order to 
avoid that the magnitude of this effect is too unprobably high, we can also push the COBE 
data to their upper (1 — 1.5)o" error bar, an increase of about (20-25)% resulting in the same 
increase in the power spectrum P{k) and a corresponding decrease of the biasing factor h (x^ 
is unaffected). We would still get acceptable numbers, though h would become rather low 
(~ 1.5), in particular ^(IO-'^^Mq), (t(1O^^M0) which are otherwise a little bit low get better. 
Note however, that even for the mean COBE normalization, the values of the large-scale 
bulk velocities in the n = 2 model for the allowed range of parameters are slightly higher 
than those in the CDM model with a cosmological constant and flat initial perturbation 
spectrum (see e.g. [0) which are in turn higher than those in the tilted CDM model with 
ns < 0.9. Thus, this problem is less severe for the double inflationary n = 2 model than for 
the latter models. 

In conclusion, we have considered here two double-inflationary models with two nonin- 
teracting scalar fields: a light scalar field with potential \mf(j)f and a heavy scalar field (ph 
with potential ^(^/J. We have analized numerically the cases n = 2 and n = 4. Trying CDM 
with an initial spectrum produced by these models, one has three free parameters, that is 
two more than the standard CDM model with a (approximately) scale-invariant initial spec- 

trum. For given /c^ (location of the "step") and given p = ^ — , normalization to COBE 
DMR constrain the remaining free parameter of the model: for n = 2, rrih ~ 3 x 10~^Mp, 
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a result already reported earlier for n = A, ^/X ~ 2 x 10 ^. In this way also, the 



energy scale of the inflationary phase is determined, scales ~ cross the horizon (for the 
first time) close to the beginning of the second inflation corresponding to an energy density 
~ ^ X 10~^^Mp for n = 2 and ~ ^ x 10~^^Mp for n = 4. As was recently emphasized in the 
study of another double-inflationary model pBf , here too the peculiar velocities obtained are 
very small for all values of the parameters and this poses a severe constraint on the model. 
The n = 4 model is shown to be excluded while the n = 2 model is marginally admissible 
for the range of parameters |^ ~ (6 — lO)h~^Mpc and 10 < p < 15. In the latter case, the 
remaining difficulty is still with low large-scale bulk velocities though it is less severe than 
in the CDM+A model or the tilted CDM model. 

If one assumes, invoking cosmic variance, that the average values of the measured peculiar 
velocities are higher than their actual r.m.s. values, and taking the upper (1 — 1.5)a limit 
of the COBE DMR data for the overall normalization, a window of parameters mentioned 
above for the n = 2 model remains compatible with observations, however with velocities 
about 30% less than the lower la measured ones. If the COBE DMR measurements go a 
little bit up as they did recently, while the measured peculiar velocities on which there are 
still rather large uncertainties godown, then these models will do better. Another potential 
difficulty which might also be cured by an overall increase of amplitude is the rather small 
total density fluctuations at galaxy scales. 
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FIGURES 



FIG. 1. Power spectrum for n = 2 and, from top to bottom on the left, p =12, 10, 15 and 28 
(full lines), compared to a scale invariant spectrum (dashed line) while all spectra are normalized 
to COBE. 

FIG. 2. Power spectrum for n = 4 and, from top to bottom on the left, p =3, 8 and 16 (full 
line), again with a scale invariant spectrum (dashed) and all spectra normalized to GOBE. 

FIG. 3. Schematic representation of the behaviour of the n = 2 model in the p — kj, 
plane with normalization according to the mean GOBE data. For p < 8, 6 < 1.5 while for 
p > 14, (t(1O^^M0) < 1.5. Above the dashed line, vqq < 215km/s. Inside the upper, resp. lower 
curve, < 7, resp. 3. 
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TABLES 



p 


8 


10 


12 


14 


15 


28 


b 


1.53 


1.78 


1.85 


1.94 


1.89 


1.63 




269 


260 


263 


263 


266 


282 


vm 


221 


215 


218 


218 


220 


231 




3.64 


2.85 


2.32 


1.99 


1.88 


1.4 




2.69 


2.1 


1.75 


1.50 


1.44 


1.26 




17.23 


8.98 


4.35 


2.08 


2.13 


4.47 



TABLE L Values of the various tests for the model with n = 2. All these number have been 
calculated with k^^ = lh~^ Mpc. 



p 


3 


5 


6 


8 


16 


b 


1.16 


1.56 


1.78 


2.2 


3.09 




298 


250 


235 


217 


208 


ve,o 


241 


206 


195 


182 


176 




5.49 


3.83 


3.22 


2.38 


1.11 


(J(10l2^o) 


4 


2.82 


2.38 


1.78 


0.87 




41.88 


27.39 


23.24 


11.39 


3.85 



TABLE II. Values of the various tests for the model with n = 4. All these number have been 
calculated with k^^ = lh~^ Mpc. 



h 


1 


1.2 


1.5 


3 


b 


1.97 


2.07 


2.21 


2.64 




260 


255 


250 


228 




215 


212 


208 


193 




2.14 


2.15 


2.16 


2.25 




1.6 


1.58 


1.58 


1.68 


x' 


2.51 


1.85 


1.97 


3.4 



TABLE III. n = 2 model for p = 13 and values of kj^^, expressed in ^ ^ Mpc, around the 
window of best parameters. 



/ 


10 


15 


20 


25 


30 


40 


50 


60 


75 




1.346 


0.773 


0.488 


0.327 


0.230 


0.126 


.0738 


.0486 


.0272 




1.24 


0.74 


0.49 


0.37 


0.24 


0.14 


0.080 


0.048 


0.025 



TABLE IV. Comparison between predicted cr^(Z), where I is the cell size expressed in /i~^Mpc 
and cr^j,^(Z), the values inferred from the Stromlo-APM redshift survey, for n = 2, p = 13 and 
k'^^ = lh~^ Mpc. For the values displayed, the test gives = 2.51. 
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